Theoretical Justification of Structural, Magnetoelectronic and Optical Properties in QFeO3 (Q = Bi, P, Sb): A First-Principles Study

One of the primary objectives of scientific research is to create state-of-the-art multiferroic (MF) materials that exhibit interconnected properties, such as piezoelectricity, magnetoelectricity, and magnetostriction, and remain functional under normal ambient temperature conditions. In this study, we employed first-principles calculations to investigate how changing pnictogen elements affect the structural, electronic, magnetic, and optical characteristics of QFeO3 (Q = Bi, P, SB). Electronic band structures reveal that BiFeO3 is a semiconductor compound; however, PFeO3 and SbFeO3 are metallic. The studied compounds are promising for spintronics, as they exhibit excellent magnetic properties. The calculated magnetic moments decreased as we replaced Bi with SB and P in BiFeO3. A red shift in the values of ε2(ω) was evident from the presented spectra as we substituted Bi with Sb and P in BiFeO3. QFeO3 (Q = Bi, P, SB) showed the maximum absorption of incident photons in the visible region. The results obtained from calculating the optical parameters suggest that these materials have a strong potential to be used in photovoltaic applications.


Introduction
Multiferroic (MF) materials have attracted the interest of researchers due to their technological applications as ferroelastic, ferroelectric, and (anti)ferromagnetic materials in a single crystalline phase that gives these materials added advantages over others.Specifically, magnetoelectric (ME) coupling is the best example of multiferroism [1][2][3].Therefore, magnetic ordering can be produced in these materials by an electric field, and electric polarization can be produced by a magnetic field [4,5].The existence of ME coupling is the primary condition for the existence of multiferroism; however, a material cannot be categorized as an MF material even if both properties are present in the material.These materials can act as potential candidates for technological applications in biomedicines, data storage devices, sensors, magnetic readers, actuators, spintronics, information processing and storage, tunneling, and others [6][7][8][9][10][11][12][13][14][15][16].Unfortunately, magnetic and ferroelectric behavior exists in few compounds above room temperature.BiFeO 3 (BFO) is the most capable multiferroic material as it shows both properties simultaneously.In recent years, numerous research articles have been published on multiferroic BiFeO 3 (bismuth ferrite) due to its extraordinary nature of having ferroelectric and (anti)ferromagnetic properties at room temperature [9,[17][18][19].Researchers have focused on single-phase BiFeO 3 because of its high antiferromagnetic Néel temperature (up to 640 K) and ferroelectric transition temperature (1100 K) [9,[17][18][19].However, feeble magnetic properties [20] and high leakage current value [21] are some of the main challenges associated with this material.The creation of secondary phases (Bi 2 Fe 4 O 9 and Bi 25 FeO 40 ) and a narrow window for thermal stability are the main hurdles in the development of unpolluted single-phased BiFeO 3 [22,23].
Recent work on MF compounds and particularly ME coupling explains ways to improve these features through the development of double perovskites [41,42], heterojunction [43,44], doping [45,46], co-doping [47,48], and modulation of different properties [49,50].In these cases, theoretical methods have been employed to address the magnetic, optical, electronic, and structural properties of the investigated compounds; however, experimental studies have explained the improvement in the magnetoelectric coupling of the investigated compounds.In this manuscript, density functional theory (DFT) calculations were used to study the structural, electronic, optical, and magnetic properties of QFeO 3 (X = Bi, P, Sb).

Materials and Methods
The main motivation to perform this study was to use the FP-LAPW+lo (full-potential linearized augmented plane wave plus local orbital) method [51,52] to obtain the solution to the Kohn-Sham equations with generalized gradient approximation (GGA) as executed in WIEN2K (version 14.0) code [53,54] in the context of DFT for perovskite materials QFeO 3 (Q = Bi, P, Sb).In these calculations, the exchange correlation potential was treated by the GGA approximation proposed by Perdew Burke and Emzerhof [55].It is a well-known fact that the optical and electronic properties of systems with strongly correlated electrons, especially those with incomplete 3d, 4f, and 5f orbitals, cannot be calculated accurately by simple GGA or local density (LDA) approximation [56].Researchers have proposed numerous theoretical methods to overcome this problem.In this manuscript, the problem of strongly correlated potentials was overcome by adding the Hubbard Hamiltonian to GGA functions (GGA+U).In these calculations, a set of U values (3.0 eV, 5.0 eV, 7.0 eV, and 9.0 eV) was employed to obtain optimized results.We chose to present results calculated with U = 5.0 eV from the mentioned list of energies.The calculated TDOS spectra for QFeO 3 (Q = Bi, P, Sb) at all U values are presented in the manuscript.The rest of the results, including magnetic properties, energy band structures, and optical properties, are presented only for U = 5.0 eV.
Unit cells were divided into two parts when working with the FP-LAPW method, i.e., the interstitial region (IR) and atomic spheres (muffin tins).The IR is a region other than the muffin tin spheres.Two completely diverse basis sets expand the wave function in these regions.The wave functions in the muffin tins and IRs are expanded with the help of the atomic-like function, as shown in Equation ( 1), and the plane wave basis, as shown in Equation (2).
Here, → k and V k are the Fermi wave vectors in the irreducible Brillouin zone (IBZ) and the scalar amplitude of the plane wave, respectively.The relativistic and semi-relativistic methods are used for core and valance electrons, respectively.Choosing the best input factors is essential in first-principles calculations.Numerous input factors were set in such a way that they yielded the most precise results and ensured excellent convergence of total energy.Suitable values of R MT were used to ensure the convergence of total energy and to prevent the leakage of charge from the core states of these compounds.The following R MT values were selected: 1.63, 1.89, 2.5, 2.37, and 2.5 for O, Fe, Bi, P, and Sb, respectively.The values of l max , G max , and R MT × k max were 10, 12, and 7.0 respectively.Here, l max , R MT , and k max were used for the maximum value of angular momentum, the smallest muffin tin radius, and the maximum value of the Fermi wave vector, respectively.Selfconsistent field calculations for the optical and electronic properties of QFeO 3 (Q = Bi, P, Sb) were performed using a k-mesh of 500 points with a modified tetrahedron approach.The rhombohedral crystalline structures of QFeO 3 (Q = Bi, P, Sb) with the space group R3c are presented in Figure 1.
Micromachines 2023, 14, x FOR PEER REVIEW of the atomic-like function, as shown in Equation ( 1), and the plane wave basis, as s in Equation (2).
Here,  ⃗ and  are the Fermi wave vectors in the irreducible Brillouin zone and the scalar amplitude of the plane wave, respectively.The relativistic and semi-re istic methods are used for core and valance electrons, respectively.Choosing the b put factors is essential in first-principles calculations.Numerous input factors were such a way that they yielded the most precise results and ensured excellent conver of total energy.Suitable values of  were used to ensure the convergence of tot ergy and to prevent the leakage of charge from the core states of these compound following  values were selected: 1.63, 1.89, 2.5, 2.37, and 2.5 for O, Fe, Bi, P, an respectively.The values of  ,  , and  ×  were 10, 12, and 7.0 respect Here,  ,  , and  were used for the maximum value of angular momentu smallest muffin tin radius, and the maximum value of the Fermi wave vector, respect Self-consistent field calculations for the optical and electronic properties of QFeO3 (Q P, Sb) were performed using a k-mesh of 500 points with a modified tetrahedro proach.The rhombohedral crystalline structures of QFeO3 (Q = Bi, P, Sb) with the group R3c are presented in Figure 1.

Electronic Band Structures
Studying energy band structure is very useful, as numerous physical propert crystalline solids are associated with energy band structure.Therefore, knowledge ergy band gaps is crucial because of dexterous applications in electromagnetic, opt tronic, and magneto-optic devices.Materials with collinear conduction band m (CBM) and valance band maxima (VBM) are known as direct band gap semicondu and others are known as indirect band gap semiconductors.The response of indirect gap materials to optical excitations is very weak compared with that of direct ban materials.In addition to photons and electrons, phonons also participate in optical t tions to preserve the momentum in indirect band gap materials, and this is the reaso researchers prefer direct band gap materials instead of indirect band gap materia applications in optoelectronic devices.
Energy band calculations were performed for QFeO3 (Q = Bi, P, Sb) along the h symmetric axis (R, Г, X, M, Г) in the irreducible Brillouin zone (IBZ) for both majorit

Electronic Band Structures
Studying energy band structure is very useful, as numerous physical properties of crystalline solids are associated with energy band structure.Therefore, knowledge of energy band gaps is crucial because of dexterous applications in electromagnetic, optoelectronic, and magneto-optic devices.Materials with collinear conduction band minima (CBM) and valance band maxima (VBM) are known as direct band gap semiconductors, and others are known as indirect band gap semiconductors.The response of indirect band gap materials to optical excitations is very weak compared with that of direct band gap materials.In addition to photons and electrons, phonons also participate in optical transitions to preserve the momentum in indirect band gap materials, and this is the reason that researchers prefer direct band gap materials instead of indirect band gap materials for applications in optoelectronic devices.
Energy band calculations were performed for QFeO 3 (Q = Bi, P, Sb) along the highly symmetric axis (R, Γ, X, M, Γ) in the irreducible Brillouin zone (IBZ) for both majority and minority spins.Calculated energy band structures with U = 5.0 eV are presented in Figure 2 for QFeO 3 (Q = Bi, P, Sb).In Figure 2, the dotted line shows the Fermi level E F , and the region below the dotted line is known as the valance band, whereas the region above the dotted line is known as the conduction band.The calculation reveals that QFeO 3 compounds (Q = Bi, P, Sb) are indirect band gap semiconductors.The calculated values of the energy bandgaps for QFeO 3 (Q = Bi, P, Sb) for various values of U are tabulated in Table 1.The different behaviors of the calculated energy band structures for the majority and minority spins showed spin anisotropy, which was responsible for the magnetic moments in QFeO 3 (Q = Bi, P, Sb).
Micromachines 2023, 14, x FOR PEER REVIEW 4 of 16 minority spins.Calculated energy band structures with U = 5.0 eV are presented in Figure 2 for QFeO3 (Q = Bi, P, Sb).In Figure 2, the dotted line shows the Fermi level  , and the region below the dotted line is known as the valance band, whereas the region above the dotted line is known as the conduction band.The calculation reveals that QFeO3 compounds (Q = Bi, P, Sb) are indirect band gap semiconductors.The calculated values of the energy bandgaps for QFeO3 (Q = Bi, P, Sb) for various values of U are tabulated in Table 1.The different behaviors of the calculated energy band structures for the majority and minority spins showed spin anisotropy, which was responsible for the magnetic moments in QFeO3 (Q = Bi, P, Sb).The density of states (total (TDOS) and partial (PDOS)) calculated by spin-polarized (SP) calculations for QFeO3 (Q = Bi, P, Sb) are discussed in this section of the manuscript.The calculated TDOS spectra for QFeO3 (Q = Bi, P, Sb) at all U values (3.0 eV, 5.0 eV, 7.0 eV, and 9.0 eV) are presented in Figure 3.A remarkable difference was present in the spectra of spin-polarized TDOS QFeO3 (Q = Bi, P, Sb) for both majority and minority spins presented in Figure 3 from −6.0 to 6.0 eV.Orbitals below −6.0 eV are known as core states, and they do not have a significant impact on the physical properties of compounds, which

Density of States
The density of states (total (TDOS) and partial (PDOS)) calculated by spin-polarized (SP) calculations for QFeO 3 (Q = Bi, P, Sb) are discussed in this section of the manuscript.The calculated TDOS spectra for QFeO 3 (Q = Bi, P, Sb) at all U values (3.0 eV, 5.0 eV, 7.0 eV, and 9.0 eV) are presented in Figure 3.A remarkable difference was present in the spectra of spin-polarized TDOS QFeO 3 (Q = Bi, P, Sb) for both majority and minority spins presented in Figure 3 from −6.0 to 6.0 eV.Orbitals below −6.0 eV are known as core states, and they do not have a significant impact on the physical properties of compounds, which means that core electrons are not used to determine magnetic, optical, thermal, and electronic properties.Likewise, orbitals present between −6.0 to −4.0 are known as semi-core electrons, and they do not have a significant impact on the physical properties of compounds either.One can note that the DOS values of the entire energy range are different for majority and minority spins, which shows that these materials possess non-zero total magnetic moments.The magnetic moment for QFeO 3 (Q = Bi, P, Sb) above the Fermi level (conduction band) is present due to the disproportion of spin in the minority spin, as the DOS in the minority spin is dominant compared with the DOS in the majority spin.
Micromachines 2023, 14, x FOR PEER REVIEW 5 of 16 means that core electrons are not used to determine magnetic, optical, thermal, and electronic properties.Likewise, orbitals present between −6.0 to −4.0 are known as semi-core electrons, and they do not have a significant impact on the physical properties of compounds either.One can note that the DOS values of the entire energy range are different for majority and minority spins, which shows that these materials possess non-zero total magnetic moments.The magnetic moment for QFeO3 (Q = Bi, P, Sb) above the Fermi level (conduction band) is present due to the disproportion of spin in the minority spin, as the DOS in the minority spin is dominant compared with the DOS in the majority spin.The spectra of the PDOS presented in Figure 4 show that the major contribution to the valance band of BiFeO3 came from O and Fe atoms, and some minor contributions from Bi atoms were also present.One can notice from the spectra of PDOS that the O-2 , Bi-6 , and Bi-6 orbitals had major contributions to the valance band, whereas the Fe-3 orbitals had some minor contributions.In the conduction band, major contributions were present due to the Fe atoms in the minority spin channel, and some minor contributions from the Bi and O-atoms were also present.One can notice from the spectra of the PDOS that Fe-3 orbitals had major contributions to the conduction band, whereas the Bi-6 and O-2 orbitals had some minor contributions.The spectra of the PDOS presented in Figure 4 show that the major contribution to the valance band of BiFeO 3 came from O and Fe atoms, and some minor contributions from Bi atoms were also present.One can notice from the spectra of PDOS that the O-2p 4 , Bi-6s 2 , and Bi-6p 3 orbitals had major contributions to the valance band, whereas the Fe-3d 6 orbitals had some minor contributions.In the conduction band, major contributions were present due to the Fe atoms in the minority spin channel, and some minor contributions from the Bi and O-atoms were also present.One can notice from the spectra of the PDOS that Fe-3d 6 orbitals had major contributions to the conduction band, whereas the Bi-6p 3  and O-2p 4 orbitals had some minor contributions.
The spectra of the PDOS presented in Figure 5 show that the major contribution to the valance band of PFeO 3 came from O-and P-atoms, and some minor contributions from Fe atoms were also present.One can notice from the spectra of the PDOS that the O-2p 4 and P-3p 3 orbitals had major contributions to the valance band, whereas Fe-3d 6 orbitals also had some minor contributions i.Major contributions in the conduction band were present due to P and Fe atoms, whereas minor contributions from O atoms were also present.One can notice from the spectra of the PDOS that in the conduction band, the P-3p 3 and Fe-3d 6 orbitals had major contributions and minor contributions from O-2p 4 .
The spectra of the PDOS presented in Figure 6 show that the major contribution in the valance band of SbFeO 3 came from O atoms, and some minor contributions from Sb and Fe atoms were also present.One can notice from the spectra of the PDOS that the O-2p 4 orbitals had major contributions, whereas the Sb-5p 3 and Fe-3d 6 orbitals also had minor contributions in the valance band.In the conduction band, major contributions were present from Fe and Sb atoms in the minority spin channel, and some minor contributions from O atoms were also present.One can notice from the spectra of the PDOS that Fe-3d 6 , Sb-5s 2 , and Sb-5p 3 orbitals had major contributions, whereas O-2p 4 orbitals had some minor contributions in the conduction band.The spectra of the PDOS presented in Figure 5 show that the major contribution to the valance band of PFeO3 came from O-and P-atoms, and some minor contributions from Fe atoms were also present.One can notice from the spectra of the PDOS that the O-2 and P-3 orbitals had major contributions to the valance band, whereas Fe-3 orbitals also had some minor contributions i.Major contributions in the conduction band were present due to P and Fe atoms, whereas minor contributions from O atoms were also present.One can notice from the spectra of the PDOS that in the conduction band, the P-3 and Fe-3 orbitals had major contributions and minor contributions from O-2 .The spectra of the PDOS presented in Figure 5 show that the major contribution to the valance band of PFeO3 came from O-and P-atoms, and some minor contributions from Fe atoms were also present.One can notice from the spectra of the PDOS that the O-2 and P-3 orbitals had major contributions to the valance band, whereas Fe-3 orbitals also had some minor contributions i.Major contributions in the conduction band were present due to P and Fe atoms, whereas minor contributions from O atoms were also present.One can notice from the spectra of the PDOS that in the conduction band, the P-3 and Fe-3 orbitals had major contributions and minor contributions from O-2 .

Magnetic Properties
Magnetic moments were also calculated in this study with the GGA+U approach.A noticeable difference was present in the spectra of PDOS and TDOS for QFeO 3 (Q = Bi, P, Sb) near the Fermi level in different spins.This difference was the origin of magnetic behavior in the compounds above.Information regarding the electronic states with a major impact on magnetic properties can be obtained from the spectra of the PDOS.The spectra of the PDOS presented in Figures 4-6 show that the magnetic moment in the valance band was mainly present due to the spin imbalance of O atoms from −4.0 to −1.0 eV; however, magnetic moment in the conduction band was mainly present due to spin imbalance of Fe atoms in the negative region between 2.0 and 4.0 eV.The spectra of PDOS for BiFeO 3 presented in Figure 5 show that magnetic moment in the valance band was mainly present due to spin imbalance of the Bi-6p 3 , Fe-3d 6 , and O-2p 4 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3d 6  and Bi-6p 3 orbitals.
Micromachines 2023, 14, x FOR PEER REVIEW 7 of 16 The spectra of the PDOS presented in Figure 6 show that the major contribution in the valance band of SbFeO3 came from O atoms, and some minor contributions from Sb and Fe atoms were also present.One can notice from the spectra of the PDOS that the O-2 orbitals had major contributions, whereas the Sb-5 and Fe-3 orbitals also had minor contributions in the valance band.In the conduction band, major contributions were present from Fe and Sb atoms in the minority spin channel, and some minor contributions from O atoms were also present.One can notice from the spectra of the PDOS that Fe-3 , Sb-5 , and Sb-5 orbitals had major contributions, whereas O-2 orbitals had some minor contributions in the conduction band.

Magnetic Properties
Magnetic moments were also calculated in this study with the GGA+U approach.A noticeable difference was present in the spectra of PDOS and TDOS for QFeO3 (Q = Bi, P, Sb) near the Fermi level in different spins.This difference was the origin of magnetic behavior in the compounds above.Information regarding the electronic states with a major impact on magnetic properties can be obtained from the spectra of the PDOS.The spectra of the PDOS presented in Figures 4-6 show that the magnetic moment in the valance band was mainly present due to the spin imbalance of O atoms from −4.0 to −1.0 eV; however, magnetic moment in the conduction band was mainly present due to spin imbalance of Fe atoms in the negative region between 2.0 and 4.0 eV.The spectra of PDOS for BiFeO3 presented in Figure 5 show that magnetic moment in the valance band was mainly present due to spin imbalance of the Bi-6 , Fe-3 , and O-2 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3 and Bi-6 orbitals.
The spectra of the PDOS for PFeO3 presented in Figure 6 show that magnetic moment in the valance band was mainly present due to the spin imbalance of the P-3 and O-2 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3 , P-3 , and O-2 orbitals.The spectra of the PDOS for SbFeO3 presented in Figure 6 show that the magnetic moment in the valance band was mainly present due to the spin imbalance of the Sb-5 and O-2 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3 , Sb-5 , and O-2 orbitals.The calculated total and partial magnetic The spectra of the PDOS for PFeO 3 presented in Figure 6 show that magnetic moment in the valance band was mainly present due to the spin imbalance of the P-3p 3 and O-2p 4 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3d 6 , P-3p 3 , and O-2p 4 orbitals.The spectra of the PDOS for SbFeO 3 presented in Figure 6 show that the magnetic moment in the valance band was mainly present due to the spin imbalance of the Sb-5p 3 and O-2p 4 orbitals.However, the magnetic moment in the conduction band was mainly present due to the spin imbalance of the Fe-3d 6 , Sb-5p 3 , and O-2p 4 orbitals.The calculated total and partial magnetic moments are presented in Table 2.We can observe from Table 2 that Fe atoms had major contributions in all these compounds.

Optical Properties
The response of the medium in the presence of external electromagnetic radiation can be explained with the help of the complex dielectric function ε(ω).The following equation can be used to calculate ε(ω): Here, ε 1 (ω) and ε 2 (ω) are known as the real and imaginary parts of the complex dielectric function.The imaginary part ε 2 (ω) of the complex dielectric function explains the absorptive behavior of the medium and is directly associated with energy band structure.The following equation can be used to calculate ε 2 (ω): Here, m, e, E k n (k), f kn , and P v nn ′ (k) are used to specify the mass of the electron, the charge of the electron, the energy of one electron, the Fermi-Dirac distribution, and the projection of momentum dipole matrix elements, respectively.These matrix elements are described between the initial and final stages of the electric field in the v direction.The calculated spectra of ε 2 (ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 7 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of ε 2 (ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The point form where peaks start to originate is known as the threshold energy of the material.The threshold energies of ε 2 (ω) calculated from Figure 7 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 3.In both spin channels, PFeO 3 showed the maximum absorption of photons in the visible and infrared region (IR), and BiFeO 3 had the lowest absorption in the visible region.One can observe from Figure 7 that the spectra of ε 2 (ω) shifted toward higher energies as P was substituted for Sb and Bi.The origin of different peaks in the spectra of ε 2 (ω) can be identified by the spectra of the PDOS.The first peak in the spectra of ε 2 (ω) occurred due to electronic transitions from the O-2p 4 orbitals to Bi-6p 3 , P-3p 3 , and Sb-5p 3  for BiFeO 3 , PFeO 3 , and SbFeO 3 , respectively.
Micromachines 2023, 14, x FOR PEER REVIEW 9 of 16 The real part  () of the complex dielectric function explains the dispersive behavior of the medium and is inversely associated with energy band structure.The calculated spectra of  () for QFeO3 (Q = Bi, P, Sb) are presented in Figure 7 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of  () for QFeO3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The most significant quantity in the spectra of  () is its static value  (0), also known as the zero-frequency limit.The static values of the dielectric constant  (0) calculated from Figure 7 for QFeO3 (Q = Bi, P, Sb) are presented in Table 4.The following relationship can be used to relate the energy band gap and the static dielectric constant.
One can observe that values of  (0) increased as Bi was substituted for Sb and P, confirming its inverse relationship with the energy band gap.There was an increasing trend in the spectra of  () after the zero-frequency limit except for PFeO3 in the majority spin channel.The spectra of  () reached their peak values, then a sharp dip was observed in the spectra, and finally, these curves entered the negative region at certain values.The frequencies at which these curves enter the negative region are known as the plasmon frequencies of the compounds.The energy ranges at which the spectra remained in the negative region showed complete attenuation of the incident beam of photons.Table 3. Threshold energies (eV) of ε 2 (ω), K(ω), I(ω), L(ω), and σ(ω) for QFeO 3 (Q = Bi, P, Sb).The following expression can be used to calculate the real part ε 1 (ω) of the complex dielectric function using the ε 2 (ω) calculated by the Kramers-Kronig transformation: The real part ε 1 (ω) of the complex dielectric function explains the dispersive behavior of the medium and is inversely associated with energy band structure.The calculated spectra of ε 1 (ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 7 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of ε 1 (ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The most significant quantity in the spectra of ε 1 (ω) is its static value ε 1 (0), also known as the zero- frequency limit.The static values of the dielectric constant ε 1 (0) calculated from Figure 7 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 4.The following relationship can be used to relate the energy band gap and the static dielectric constant.
One can observe that values of ε 1 (0) increased as Bi was substituted for Sb and P, confirming its inverse relationship with the energy band gap.There was an increasing trend in the spectra of ε 1 (ω) after the zero-frequency limit except for PFeO 3 in the majority spin channel.The spectra of ε 1 (ω) reached their peak values, then a sharp dip was observed in the spectra, and finally, these curves entered the negative region at certain values.The frequencies at which these curves enter the negative region are known as the plasmon frequencies of the compounds.The energy ranges at which the spectra remained in the negative region showed complete attenuation of the incident beam of photons.Table 4. Static values of ε 1 (0), n(0), and R(0) r QFeO 3 (Q = Bi, P, Sb).Knowledge of ε 2 (ω) and ε 1 (ω) can be used to calculate the rest of the optical parame- ters such as the extinction coefficient K(ω), refractive index n(ω), reflectivity coefficient R(ω), absorption coefficient I(ω), energy loss function L(ω), and the real part of optical conductivity σ(ω).The following equation can be used to calculate K(ω) by utilizing the knowledge of ε 2 (ω) and ε 1 (ω): The extinction coefficient K(ω) showed features similar to those of ε 2 (ω), which means that it also tells us about the absorptive behavior of the compound.The calculated spectra of K(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 8 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of K(ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The point form where the peaks originate is known as the threshold energy of the material.The threshold energies of K(ω) calculated from Figure 8 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 3.In both spin channels, PFeO 3 showed the maximum absorption of photons in the visible and infrared region (IR), and BiFeO 3 had the lowest absorption in the visible region.One can observe from Figure 9 that the spectra of K(ω) shifted toward higher energies as P was substituted for Sb and Bi.The following expression explains the reflectivity coefficient () in terms of () and () [58]: The calculated values of the complex dielectric function () can also be used to explain the behavior of () with the help of the following expression [58]: The calculated spectra of () for QFeO3 (Q = Bi, P, Sb) are presented in Figure 9 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of () for QFeO3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The most significant quantity in the spectra of () is its static value (0) , also known as the zero-frequency limit.The static values of (0) calculated from Figure 9 for QFeO3 (Q = Bi, P, Sb) are presented in Table 4.The calculated results showed that the reflectivity () of PFeO3 was greater than that of BiFeO3 and SbFeO3.From spectra of (), it can be concluded that PFeO3 can be used in applications in which a maximum of 60% reflection of incident photons is needed.BiFeO3 and SbFeO3 are poorly reflecting materials in the entire energy range.An interesting behavior of reflectivity can clearly be seen in Figure 9, in which the peak value of reflectivity occurs in the region where  () becomes negative because the reflectivity of the compound and metallicity are directly associated with each other.The region where  () becomes negative reveals the metallic nature of compounds [59].The reflectivity for BiFeO3, PFeO3, and SbFeO3 started from 12.7%, 28.8%, and 16.6%, respectively, for the majority spin channel.However, reflectivity in the minority spin channels for BiFeO3, PFeO3, and SbFeO3 started from 15.6%, 24.8%, and 18.7%, respectively.The reflectivity attained its maximum values of approximately 42%, The loss of rapidly moving electrons through the compounds can be determined by the energy loss function ().Valuable information about the interaction of electronic systems with incident photons can be obtained from the spectra of ().The calculated spectra of () for QFeO3 (Q = Bi, P, Sb) are presented in Figure 9 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of () for QFeO3 (Q = Bi, P, Sb) are different in the minority and majority spin channels.The point form where peaks originate is known as the threshold energy of the material.The threshold energies of () calculated from Figure 9 for QFeO3 (Q = Bi, P, Sb) are presented in Table 3.The energy loss functions () and  () are inversely associated with each other, which means that () is active where  () is silent and vice versa.The continuation of electrical transportation at higher energies of photons is known as optical conductivity.Optical conductivity () is a very useful tool for the study of conducting materials.The following equation can be used to calculate (): The following equation can be used to calculate n(ω) by utilizing the knowledge of ε 2 (ω) and ε 1 (ω): The refractive index n(ω) is an important parameter to check the suitability of an optical material for its technological applications in detectors, waveguides, photonic crystals, and solar cells.The refractive index n(ω) shows features similar to those of ε 1 (ω), which means that it also tells us about the dispersive behavior of the compound.The calculated spectra of n(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 8 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of n(ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The most significant quantity in the spectra of n(ω) is its static value n(0), also known as the zero-frequency limit.The static values of n(0) calculated from Figure 8 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 4.The following relationship can be used to relate static values of the dielectric constant with the static values of the refractive index: One can observe that values of n(0) increased as Bi was substituted with Sb and P.There was an increasing trend in the spectra of n(ω) after the zero-frequency limit except for PFeO 3 in the majority spin channel.The spectra of n(ω) reached their peak values, then a sharp dip was observed in the spectra, and finally, these curves became less than unity (n = c/v) at certain energies.This shows an unnatural phenomenon in which the speed of light in vacuum (c) is less than the group velocities of the incident photons (V g ).This means that the medium shows non-linear behavior instead of linear behavior, or we can say that the medium becomes superluminal [57].
The following expression explains the reflectivity coefficient R(ω) in terms of n(ω) and K(ω) [58]: The calculated values of the complex dielectric function ε(ω) can also be used to explain the behavior of R(ω) with the help of the following expression [58]: The calculated spectra of R(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 9 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of R(ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The most significant quantity in the spectra of R(ω) is its static value R(0), also known as the zero-frequency limit.The static values of R(0) calculated from Figure 9 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 4.The calculated results showed that the reflectivity R(ω) of PFeO 3 was greater than that of BiFeO 3 and SbFeO 3 .From spectra of R(ω), it can be concluded that PFeO 3 can be used in applications in which a maximum of 60% reflection of incident photons is needed.BiFeO 3 and SbFeO 3 are poorly reflecting materials in the entire energy range.An interesting behavior of reflectivity can clearly be seen in Figure 9, in which the peak value of reflectivity occurs in the region where ε 1 (ω) becomes negative because the reflectivity of the compound and metallicity are directly associated with each other.The region where ε 1 (ω) becomes negative reveals the metallic nature of compounds [59].The reflectivity for BiFeO 3 , PFeO 3 , and SbFeO 3 started from 12.7%, 28.8%, and 16.6%, respectively, for the majority spin channel.However, reflectivity in the minority spin channels for BiFeO 3 , PFeO 3 , and SbFeO 3 started from 15.6%, 24.8%, and 18.7%, respectively.The reflectivity attained its maximum values of approximately 42%, 58%, and 45% for BiFeO 3 , PFeO 3 , and SbFeO 3 , respectively, for the majority spin channel.However, the maximum reflectivity in the minority spin channels was approximately 48%, 40%, and 46% for BiFeO 3 , PFeO 3 , and SbFeO 3 , respectively.
The following equation can be used to calculate the energy loss function L(ω) by utilizing the information of the dielectric function ε(ω) [58]: The loss of rapidly moving electrons through the compounds can be determined by the energy loss function L(ω).Valuable information about the interaction of electronic systems with incident photons can be obtained from the spectra of L(ω).The calculated spectra of L(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 9 for both the minority and majority spins.The heights and shapes of the peaks in the spectra of L(ω) for QFeO 3 (Q = Bi, P, Sb) are different in the minority and majority spin channels.The point form where peaks originate is known as the threshold energy of the material.The threshold energies of L(ω) calculated from Figure 9 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 3.The energy loss functions L(ω) and ε 2 (ω) are inversely associated with each other, which means that L(ω) is active where ε 2 (ω) is silent and vice versa.
The continuation of electrical transportation at higher energies of photons is known as optical conductivity.Optical conductivity σ(ω) is a very useful tool for the study of conducting materials.The following equation can be used to calculate σ(ω): Here, σ 1 (ω) and σ 2 (ω) are known as real and imaginary parts of optical conductivity.The following equation can be used to calculate optical conductivity σ(ω) by utilizing the information of the dielectric function ε(ω) [58]: Transitions between CB and VB are used to calculate the optical conductivity of the aforesaid compounds.The calculated spectra of σ(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 10 for both the minority and majority spins.The heights and shapes of peaks in the spectra of σ(ω) for QFeO 3 (Q = Bi, P, Sb) are different in the minority and majority spin channels.The point form where peaks originate is known as the threshold energy of the material.The threshold energies of σ(ω) calculated from Figure 10 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 3.
The following equation can be used to calculate I(ω) by utilizing the knowledge of ε 2 (ω) and ε 1 (ω) [58]: The absorption coefficient I(ω) can be used to obtain information regarding the penetration length of the photon in the material or the length traveled by the photon inside the material before its complete absorption while having an energy greater than the energy band gap.It also tells us about the absorptive behavior of the compound.The calculated spectra of I(ω) for QFeO 3 (Q = Bi, P, Sb) are presented in Figure 10 for both the minority and majority spins.The heights and shapes of peaks in the spectra of I(ω) for QFeO 3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The point form where the peaks originate is known as the threshold energy of the material.The threshold energies of I(ω) calculated from Figure 10 for QFeO 3 (Q = Bi, P, Sb) are presented in Table 3.The spectra of I(ω) showed increasing behavior with increasing energies of photons.
minority and majority spins.The heights and shapes of peaks in the spectra of () for QFeO3 (Q = Bi, P, Sb) were different in the minority and majority spin channels.The point form where the peaks originate is known as the threshold energy of the material.The threshold energies of () calculated from Figure 10 for QFeO3 (Q = Bi, P, Sb) are presented in Table 3.The spectra of () showed increasing behavior with increasing energies of photons.

Conclusions
Scientists seek to fabricate novel multiferroic (MF) materials that exhibit interconnected properties such as piezoelectricity, magnetoelectricity, and magnetostriction and can operate effectively at room temperature.The effect of substituting pnictogen elements on optoelectronic and magnetic properties in BiFeO3 was investigated using first-principles-based DFT calculations.We noted that electronic states moved toward the Fermi level as we substituted Bi with P and Sb.The electronic band structures revealed that PFeO3/SbFeO3 are metallic compounds, whereas BiFeO3 is a semiconductor.The studied compounds are promising for spintronics, as they exhibit excellent magnetic properties.The calculated magnetic moments decreased as we replaced Bi with SB and P in BiFeO3.A red shift in the values of  2 () was evident from the presented spectra as we substituted Bi with Sb and P in BiFeO3.QFeO3 (Q = Bi, P, SB) showed the maximum absorption of incident photons in the visible region.The results obtained from calculating optical

Conclusions
Scientists seek to fabricate novel multiferroic (MF) materials that exhibit interconnected properties such as piezoelectricity, magnetoelectricity, and magnetostriction and can operate effectively at room temperature.The effect of substituting pnictogen elements on optoelectronic and magnetic properties in BiFeO 3 was investigated using first-principlesbased DFT calculations.We noted that electronic states moved toward the Fermi level as we substituted Bi with P and Sb.The electronic band structures revealed that PFeO 3 /SbFeO 3 are metallic compounds, whereas BiFeO 3 is a semiconductor.The studied compounds are promising for spintronics, as they exhibit excellent magnetic properties.The calculated magnetic moments decreased as we replaced Bi with SB and P in BiFeO 3 .A red shift in the values of ε 2 (ω) was evident from the presented spectra as we substituted Bi with Sb and P in BiFeO 3 .QFeO 3 (Q = Bi, P, SB) showed the maximum absorption of incident photons in the visible region.The results obtained from calculating optical parameters suggest that these materials have a strong potential to be used in photovoltaic applications.

Figure 5 .
Figure 5. Partial density of states for PFeO3.

Figure 4 .
Figure 4. Partial density of states for BiFeO3.

Figure 5 .
Figure 5. Partial density of states for PFeO3.Figure 5. Partial density of states for PFeO 3 .

Figure 5 .
Figure 5. Partial density of states for PFeO3.Figure 5. Partial density of states for PFeO 3 .

Figure 6 .
Figure 6.Partial density of states for SbFeO3.